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ee A Further Study of 
= the Multichoinne! Waximun Sntropy Spectral Analysis 


!, Introduction 


In the field of signal detection and estimation, it is always 


| necessary to face the probleme with systems whose states are des- 
q eribed by several time variables. tor example, in the geophysical 
research, data are gathered by arrays of sensors (seismographs, 
hydrophones, antennas, meteorological stations). The sensors used 
in these array systems will-not always be along a line or ina 
regular grid but will often have some artitrary spatial distribu- 
tion. Riso, in an array system, one or more sensor: Can Ve connected 
4 to a channel; each of these chanriels provides one of the signal to be 
processed. Instead of processing the signals received or: each channel 
ceparately, it is possible to process all the received signals 
collectively by means of a multicharnel filter that take: into account 
the interrelations among tne individual sigmals. "Therefore, mul ti- 
channei time series nrovessing is quite important because of its 
sufficient informations about the object. To solve this problem, 
spectral estimation is an effective tool to describe the -isnad 
properties. Various attempts as - fal, | s | - ack of mabe tchannel 
upectral analysis have been made in the literature. Méeny conventionads 
methods have been extended from one channel to multicharnel. One 
rf technique which has been studied extensively due to its Nigh resolu- 
tion characteristics is the maximum entropy metnod (MEM). Jones [a]. 
Wuttal i4|, Strand 1 | aud Morf et al. [2 | have independent!y cbitained 


the multichsrnei generalization of the scalar Burg's miximuni entropy’ 


¥ method, In this report, we use the generalivation method of Strana 


” 


was 


as a means of study. Detuiled computer results are presented for 

the multichannel maximum entropy spectral analysis of shurt sinusoids. 
Although for the real signals the snectral line shifting and 
spontaneous line splitting are observed, we have demonstrated that 
such phenomena are almost completely eliminated in the complex 

signal case. This result is practically significant as it indicates 
that accurate spectral estimation is possihie for the cotuplex signal 
case. The multichannel maximum entropy spectral «nalysi: can thus 
provide a good phase estimate from which time delays amon; spatially 
separated sensors can be properly determined. 

This report is organized as follows. In section IT, we briefly 
review the algorithm of Strand. In section ITI, we show the results 
of the spectral estimation and the phase difference between channels. 
In section IV, we present the comparative result of the real signals 
and the complex signals in the viewpoint of line shifting and solit- 
ting properties. And in tne Appendix, we list a detaited computer 
program in Fortran which we used to obtain the complex multichannel 


maximum entropy spectral estimation. 


II. Theoretical Background 
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be a complex zero-mean, wide-sense stationary p-charmel time 


series. The N-element forward filter (or forward filver of length N) 


has the form. 


NL 
y(t) + y Fe y(t-kat) = e,,(t) 
wel 


Similarly, we consider the N-element backward filter 


* 
y(t) + Byyy (trket) = by (t) 
Ke 


In (2) and (3), the forward and backward filter coefficients, Fy and 

Bury respectively, are complex pxp matrices. The asteris« * denotes 
the (Hermitian) conjugate transpose. For the optimum rerward or 

backward filters (2) or (3), we impose the condition that the 
expected mean-square value of ent) or by(t) should be a minimum. 

For (2) this involves determining Fane Fone cee + “ype Such that 
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where I is the pxp identity matrix and the square hlock s'.omatrices 
in this system are jefined by 

Re =F ( y(t) ye(t- kat}, k=O,1,..68 (6) 
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We call the coefficient matrix in (5) the R-matrix, or simply R. 
The forward power matrix P,, for the resulting optimum filter (that 


is, for F., satisfying (5)) is 

Py = E felt) eye (t)} i 
By multiplying both sides of (5) by F.* we have, 
PN 


= Fy* R PY 


where we derine Fy by 
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with comma indicating matrix partition. 


Similarly, ror the uptimum backward filter, we can derive a system: 
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we call the coefficient matrix in (10) the R’ matrix, or simply R’. 


% The backwaru power matrix Di for the resulting ovtimum filter is 
¢ * } 
BS Be Das 5 
N E | D ay ) oy ( yy (11) 
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Purg has derived the following resuits concerning the recursive 


solution of (5) and (10), the forward and vackward filters may be 


postulated in the form 
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where Bes 17 | ene -- Ps = is the block reverse of Fi, ,. 


wid the pxp matrices C. and Ce are called forward and backward reflection covcfFictents, 
N ™ 


respectively. Lt can be sown that P,, and P! may be computed by using the modern 


N N 
Levinson alweritiom: 3 
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For the positive definite Pip O (and Py > 0), it is *nown that 
the maximum entropy (autoregressive) power Ree Has density, S(f), 
a pxp matrix, can then be calculated from tne forward filter F,; by 


S(t ye sale are dP, [poe f )| (17) 


where F(z)= 1+i,..4+ «+. 4F),.2 5 2 is the complex scalar defined by 


2= exp | -2nipal and £ is the frequency. 

It can be show (as suggested in (9)) that Foy PyesseF) and 
(rity Paty nee Py ) are positive definite if and only if R is positive 
detinite. Instead of trying to solve (5) using approximations to the 


Rie Strand Ce derive/a method cf estimating the reflection coefficients 


iy 


Cc and Ci by minimizing a certain weighted sum of squares cf 


residuals for the vorward and Lackward filters, after wnich all 


forward und backward filter coefficients Pyjy and Bie] are 
the modern Levinson algorithm. It can be shown uhat 
BC 


N Y-1 


where 


and Wn are positive scalar weights (usually all taken as 


oltained by 


+ Pi, Para Nh = -2G (18) 
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practice). Nd is the number of observation of © realization 


y(t), N is the number of tuples of data, 
Also au and 2a can be obtained from 
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rit= @ N-d - x v4 
ss ae i Copa mrt 


m=1,2,...Nd-N = 


The whole computing procedure is iljustrated in Fig. 1. 
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4 1fI. Numerical Results of the Multichannel WPM Snectrum and Estimation 
of Phase Ditferences Between Channels 
9 To ensure the correciuness of the progrum, the came data which { 
7 | 
r Strand used have been tested: 
: 
Nd (number of data points sampled) = 12% 
i 
At = 1 second. 
PE p (number of channels) = 2 
= r a 5 ap 
t V+ = 1] ’ SL aly ehetiecy 12°r. | 
: [Yee] 
Res Yip > evs (2.9) Bede DP OS. rane Boe GG 
Aq FS Par tu. (Ran ro - O35) 
| Tw Yui in ( tf ) ( ) (21) 
t Re Yo, = Cos (2 bh boot it: i: es ORR SCN ARE 085) 
Poe Yoq F Sate (2b f toe dee ae (Rant GG) 
where Rant is a pseudo-ranuom number uniformly dictributed vetween 
© and 1, Yd is the phase difference between channels and a is the 
input frequency. Jf we derote the snectruin by | 
se : ; j 
had 4" > f a ‘ 5 
; S(f) = | Pia) pot) (22) | 
where £ is the frequency in Ha, then the suvered coherence is given 
z by ig 2 
=z 2 er af tl : 
a coh’(f) = san a a seein (23) 
& 4 bf 9,41) 
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j A rough indication of the expected performance is given by Akaike's 
ea 
| FPE criterion, given by | 
Ned rpy yf 
; FPE(N) = det PL ( -* Smee ee 
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Figs. 2a and 2b show 344) and Soo(f) for the real sinusoid inputs 
with f= 0.0625 Hz, Va = lorad. and NPEF (number of terms: in the 


prediction error tilter) = 7. The results for (f) are snown (in 


4g 
amplitude and phase) in Figs. 2c and 2d. From Fig. “d, we can observe 
that the pnase of Sy. (1) appears to equal -1 radian wt f-'.0625 Hz. 


In Fig. 2e, the squared coherence is show. Fig. 2f£ und Fig. 2e also 


indicate the maximum and minimum eigenvalues of spectra. igs. 3a- 


Ww 
ae 


is the repeated procedure of Figs. Za-2g but with complex input sienais. 


4 


The comparison shows that the complex multichannel Vil vives better 
resolution and accuracy than the real case. 

It is well known that the phase difference between channels is 
a very important parameter. Correct phase difference can give very 
zood estimation of time delay so as to detect tne location of the 
target. Different sets of data have been tested using the comulex 
multichannel MEM program with phase difference A gs U.5 vrad,-1l.rad, 
z. rad. (f,= 0,0625 Ha), Vd= -0,5 rad (£,=0.2 Hz). The results of 
1-1 spectrum $4,() and phase of 1-2 spectrum, arg. Sof) are shown 
in Figs. 4(a,d), 5(a,bv), 6(a,b) and 7(a,b) respectively. 


It is noted that the results are satisfactory. 


IV. The Splitting and Shifting Problem in Multichannel MEM Power 
spectru 
Irn, 1974, Chen and Stegen (7| ran some computer experiments. 
They noted that the estimated spectral maximum was found to shil't in 
relation to the true maximum. The frequency shifts depended upon 
the initial phase and the length of the sinusoid. Pougere et, al. [6 | 


showed that tne problem arose from the assumption of Burg, that 


previously determined reflection ccefficient are held fixed in order 
to determine the present reflection coefficients. Also, ‘ougere [=| 
recently indicated that these properties oeecurres to be single 
channel complex signals and the multichannel real sirnals. In 

order to study this problem we have tested a sigiual whie', exemplified 
the problem for the single-channel signals but in both real and 
complex cases. ‘his is a single sine wave of frequency 1 Hz sampled 


every 0.05 seconds, with additive noise: 
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Re f(r) ~ cos( (n-1) 7 /10 + , ) Ag, (n} 

La fin) = sin( (n-1) 7 /10 + Y) + Ag, (n} 

Re r..(n) =ecos ( (n-1) 7 /10 + ee Az.,(a) 

Im (nm) =sin ( (n-1) 7 /10 + (2) + Az (n) 
where Wy and i, are initial phase, in each channel and fs, (n) and 
Eo (n) are independent realizations of random numbers in the 


range +0.5 and A is the noise amplitude. Also, we used hd and Nv@h 


to represent the number of data points and the number of terms i: 


= the prediction error filter respectively. 

7 we have noted that the number of terms of filter (NFFF) affect 
the quality of the spectrum as indicated hy Fig. a with NPEF=3,9,21. 
F The same data are applied to the real and complex multichannel 
7 cases as shown in Fig. 8b, 8c respectively. We can see the real 

| multichannel case is worse even than the single channel. case (for 

J NPEF=21), hut the compiex signal performs goud. Figs. 9-11 show 

! the eftects of the number of data sample (lid) to the location of 
spectral peak of real signals. Figs. 9a, 10a, lia show the loca- 
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tion of spectral peak is farther from 1 liz as the cycles are not 


nn —_—_—- war renee tahitee line cated. AE eat all 
ae yes aes - : acd suet ms 8 a ROR E Lats . SoA ees ta Sie - omhi t. ae 


4 
Ef 
* 
‘ 


ae. ‘ ecta 042 Op again et ota. ig oe Rebbe ie it ea si iE ae a mae He ee 
a 


wa 


completed. When the number of data points increares, the trequency 


snift decreases ay shown in Fig. Tb. AS expected there th ee node et. 
Nd#21, 41, ..., every Lull cycle, where there is excelrent agrermernt 


between true and apparent frequencies. lHiowever, “Mipu. dea, ja show 


two more cases but with different initial phase angsies. © lecations 
of spectral peak for various number of data are also wnhoe: in 
Figs. 15 and 16 respectively (for 245, A ao Seals AE leh 

Y, =-135). Again we huve noted the came nrorevty while J. atrected 
by the number oi data points sampleg. Comparing tie Loca’ ion cf 


the peak of Nd=21 of Figs. Ga, Téa and 138, and the verticnh axis ul 


Figs. 14, 15, and 16, we have noted the peak of tre real mustichanie? 


case is also atfected Ly ibe initial phase. “he posak with 57, 135" 
: . : ) Sao Sijedh ce, tA 
phases shifts farther than with oO’, 180° case. This is the same 


condition as in the single channel case. 

A very interesting property is that the location of stectral 
peaks does not shift away from the true value Yor varict: vuriations 
(number of data sampled, mumber of terms in the prediction errev 
rilter, value of initial phase). From Figs. 9b-13b und Firs. 14-16, 
the complex multichannel case always gives very satisfactory results. 
Therefore, we can conclude that the complex MEM nower spectrum 
estimation is not affected hy the factors indicated above because 


or its sufficient informations. 
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